Dewetting of a solid monolayer 
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We report on the dewetting of a monolayer on a solid substrate, where mass transport occurs 
via surface diffusion. For a wide range of parameters, a labyrinthine pattern of bilayer islands is 
formed. An irreversible regime and a thermodynamic regime are identified. In both regimes, the 
velocity of a dewetting front, the wavelength of the bilayer island pattern, and the rate of nucleation 
of dewetted zones are obtained. We also point out the existence of a scaling behavior, which is 
analyzed by means of a geometrical model. 



Liquid films, once spread on a substrate, may break- 
up into droplets to lower the surface energy. Such a pro- 
cess is called dewetting. As for liquids, thin solid films 
may break-up into droplets. However two main differ- 
ences may be pointed out. Firstly, solids exhibit strong 
surface anisotropy whereas liquids are usually isotropic. 
Secondly, mass transport mainly occurs via surface dif- 
fusion on solids at small scales, while it is mediated by 
hydrodynamics in liquids. 

Dewetting of solid layers with sub-micron thicknesses 
was observed in recent experimental studies 
Spontaneous breakup of the film into dots can be an- 
alyzed within the frame of continuum models including 
an effective wetting potential, with surface energy [4], and 
elastic effects [5]. Moreover, the nonlinear dynamics of the 
edges of these layers may also lead to the periodic 

formation of holes behind the dewetting rim. 

For even thinner films, such as Inm thick Ag on Si[8], 
one expects the discreteness of the underlying crystalline 
lattice to come to the fore. In order to investigate these 
effects, we study the dewetting of the thinnest possible 
layer: a monolayer. In order to focus on the basic pro- 
cesses, we discard effects related to substrate roughness, 
elastic interactions, or alloying. We focus on the case 
where dewetting occurs via the nucleation of holes, sub- 
sequently invading the whole film. This occurs in a well 
defined temperature window: if the temperature were 
too low the surface would be frozen; if it were too high 
-above the roughening transition- the film would be un- 
stable and would break up into a microscopically disor- 
dered pattern. 

We show that monolayer dewetting proceeds differ- 
ently from thicker layers dewetting. As shown on Fig.l, 
monolayers initially lead to a labyrinthine pattern of bi- 
layer islands, which then slowly thicken into 3-layer, and 
then 4-layers islands, etc. Two different regimes for 
monolayer dewetting, henceforth denoted as regimes I 
and II, are analyzed. While both regimes exhibit the 
same temporal scaling behavior, their microscopic dy- 
namics is qualitatively different. 

We employ Kinetic Monte Carlo (KMC) simulations 
in order to mimic experiments with a minimum number 



of ingredients. We use a Solid on Solid model on a 2D 
square lattice, with lattice unit a, and local height > 0. 
The substrate surface, at = 0, is perfectly flat and 
frozen. Epilayer atoms hop to nearest neighbor sites with 
the rates Vn when h = 1, and Un when h > 1, with 



%J/T^Es/T. 



T,J/T 



(1) 



where uq is a constant frequency, T is the temperature (in 
units with = 1), n is the number of in-plane nearest 
neighbors, J is the bond energy, and Es is the adsorbate- 
substrate interface energy. When Eg/ J is large, the ad- 
sorbate minimizes the total energy by creating high is- 
lands with steep sides [19]. When Eg 0, the adsorbate 
wets the surface. We use periodic boundary conditions, 
and lattice sizes ranging from 200 x 300 to 10^ x 10^. 
The initial condition is either: (ICl) a complete mono- 
layer with h = 1 everywhere, or (IC2) a monolayer with 
uncovered substrate {h = 0) in a straight stripe (of width 
10a). In all cases, the dewetting front invades the whole 
substrate, leaving a labyrinth of bilayer islands with a 
width A in a denuded zone behind. In IC2, with small 
enough Es and T, the dewetting front starts from the 
existing monatomic step and invades the monolayer at 
constant velocity V, while in the other cases, holes are 
nucleated and closed fronts expand in the film. 

Regime I is observed for low T/ J, and moderate Es / J. 
On a mono-atomic step which separates two bilayer is- 
lands, the density of kinks is then very low so that most 
of the time, there is no kink. Step motion then starts 
with the detachment of an atom with n = 3, at the rate 
ra, followed by a fast kink-zipping process, during which 
atoms at kink sites detach at a rate r2 ^ rs. [20] Since 
the step undergoes a total displacement a at each zip- 
ping, the velocity is equal to the total detachment rate 
from a step segment of length A 



V ~ Arq. 



(2) 



Since T/ J is small, islands grow in a irreversible fash- 
ion once a dimer is formed on the second layer. The 
rate of formation of dimers confined between two islands 
in a area ~ A^ is ~ D\^n{^ where ni is the adatom 
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FIG. 1: KMC simulations of monolayer dewetting in regimes 
I and 11. Lighter means higher, (a): IC2, (b): ICl, early 
dynamics; (c): ICl late dynamics. Es/ J — 0.75, J/T = 5 in 
(Ia,b,c). Es/J = 0.1, J/T = 2.5 in (Ila), and Es/J = 0.12, 
J/T = 0.5 in fllb.cV 
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FIG. 2: Dewetting front velocity V and width A of bilayer is- 
lands in regime I, with IC2. Dashed lines correspond to Eq.(|2|) 
with prefactor 0.15, and solid lines to Eq.JS) with prefactor 
5. Symbols are simulations. A: for J/T — 10 •, Es/T = 6 o. 
V/X: for J/T = 10 filled A, Es/T = 6 □. 

density, and D = vqo? is the diffusion constant on the 
second layer. Such a mean field expression of uo is valid 
up to logarithmic corrections [9]. Adatoms are produced 
in bursts during zipping processes, which last during a 
time tz ~ A/ar2, and ni ~ r2/D. The probability of 
dimer formation during the zipping process is therefore: 
Pz = ootz ^ X^r2/aD. For self-consistency, islands must 
be formed after the step has moved over a distance A, so 
that {X/a)Pz ~ 1, leading to 

A/a ~ r-^'^ ~ ^JnT-Es/iT 

As seen in Fig 121 Eqs.([2|) and (j3|) are in agreement with 
KMC simulations. 



In regime II, temperature is higher, and mass exchange 
between mobile atoms and steps is fast, so that mass 
transport is limited by diffusion. Since the typical dis- 
tance between the edge of the monolayer and the edge of 
the 2D islands on top of it is A, one has 

^-^(c«-4^,^) (4) 

where Ceg is the concentration of adatoms on top of the 
monolayer in equilibrium with the atomic step separating 
h = i — 1 and h = i. From detailed balance, c^ej ~ 
e(-2^+^^)/^, and cg^ ^ e-^^/^. 

The separation A is determined by the standard nu- 
cleation theory, valid at high temperatures. The Gibbs 
energy of a circular nucleus of radius R is the same for 
islands or holes: 

G = 27rRj - TiR^Es/a^ . (5) 

where 7 is the line tension. Since OrG < for > 
Rc = o?^ / Es^ nuclei with R> Rc grow irreversibly. The 
island nucleation process in the vicinity of the dewetting 
front starts with a concentration c = c^J in a typical 
area A^, and ends with c = c^q in equilibrium with the 
newly formed island of radius Rc. Mass conservation then 
imposes that X^ci'q = A^Ceg^ +7ri?^/a^, so that: 

(2) 

This relation is also valid for holes (we start from Ceg , 
and end with Ceg\ and a hole of radius Rc). The island- 
hole nucleation symmetry is an important difference with 
the low temperature case. In simulations on Fig.l Ila and 
1 lib, we indeed observe that both 2D islands and holes 
form in the vicinity of the dewetting front. 

The front velocity V and the island width A obtained 
by simulations with IC2 agree well with theoretical ex- 
pectation Eq. dl]) and (|6|), as shown in Fig. 3, at small 
Es. For large Es^ holes will be nucleated in the bulk 
of a monolayer and evaluation of V and A is hindered. 
Therefore, we have also performed simulations where the 
motion of atoms with n = 4 is forbidden. This prevents 
the formation of holes in the film. Figl3] shows that sim- 
ulation results with and without holes are in agreement 
with Eqs.(|4j) and (|6]). For the calculation of A , w e have 
used the expression of 7 from the Ising model [lO[. 

Having a fairly complete description of the motion of 
the dewetting front in regimes I and II, we shall now 
proceed with the analysis of nucleation and growth of 
Dewetted Zones (DZ) by means of a simple geometric 
model. In this model, DZs are nucleated with the rate 
J and expand isotropically in the monolayer. The pe- 
riphery of a DZ is initially a circle expanding with the 
constant velocity V, which was determined above. While 
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FIG. 3: Regime II: (a) VX is a linear function to e^^^^ — 1, 
at J/T = 2.86, IC2. Squares are normal simulations with 
Es small. Circles are simulations where hole nucleation is 
forbidden, (b) A for T/J = 0.35, and IC2, varying Es/J (o); 
T/J = 0.4, and ICl, varying Es/J (o); Es/T = 0.4, and IC2, 
varying J/T (A). Lines are from Eq.jG]) with prefactor 1. 

they grow, DZs will meet. The sum C of the perimeters 
of the DZs if they did not see each other (i.e. the sum of 
the lengths of the circles) obeys 

dtjC = 27tVN, (7) 

where is the number of DZs. Moreover, new DZs can 
only form outside pre-existing DZs, so that 

dtN = J{A-ADz), (8) 

where A is the total area, and Adz is the total area of 
the DZs. Because of the overlap of the growing DZs, 
the length Cdz of the frontier between the DZs and the 
monolayer is smaller than £. In a mean field picture, we 
assume that Cdz ^ >C(1 — Adz/A)^ so that 

dtAoz = CdzV ^ C{1 - Adz/A)V . (9) 

Instead of analyzing the full solution of ()7l8l9p , we shall 
focus on the behavior of the uncoverage = {1/2) Adz / A 
(the fraction of the substrate surface which is not cov- 
ered). At short times, one finds: 

= 7TV^Jt^/6, (10) 

and at long times 6 1/2. The saturation time tgat^ 
which separates these two regimes, may be defined as 
the time where \dttO\ is maximum. We find: 

t^,, = (11) 

where a ^ 0.85. The time tsat may be interpreted as 
the time needed for the dewetting process to invade the 
whole monolayer. 

As shown on FigHl the t^-scaling behavior of 0{t) is 
in agreement with the simulations, except for very short 
times, where originates from the thermal roughening 




FIG. 4: Uncoverage ^ as a function of time (in Monte Carlo 
Steps per Site). Simulations with: J/T — 5, Es/J = 1, □, 
and J/T = 2.5, Es/J = 0.15 o. Dashed lines are ^ t^, and 
solid lines are the solution of the geometrical model. 

of the monolayer. In the following, we shall see how tgat 
from the simulations compares with Eq. (pT]) . To do so, 
we need to derive the expression of J. 

In regime I, we resort to rate equations to describe 
the nucleation process. A hole of size N is denoted as 
a N-hole. The initial 1-hole concentration pi essentially 
results from adatom-pair formation, at rate a~^r4, and 
annihilation, at the rate Dnipi^ so that 

Let us denote cfn^n±i the rate at which a A/'-hole trans- 
forms into a (A^ ± l)-hole. The rates of size-decreasing 
transitions are related to the attachment of a mobile 
atom: (12,1 ~ 0-3,2 ~ (JA.d, ^ Dni. Other rates re- 
sults from the detachment of a 3-neighbors atom: (Ji^2 ~ 
(72,3 ~ cr4,5 ~ ^3, or 2-neighbors atoms: (73,4 ~ (75,6 ~ ^2- 
We only consider compact square 4-holes, because they 
are the most stable 4-holes [11]. As a consequence, we 
also only consider L-shaped 3-holes, which are the ones 
leading to square 4-holes. Assuming a steady-state, and 
following the standard Becker-Doring analysis, we obtain 
the hole formation rate: 

j=pif— +^+^2i^+....y\ (13) 

V^12 cri2Cr23 cri2Cr23Cr34 / 

which leads to 

<721 

when J < Es < 2 J. This regime correspond to a critical 
size of 4. Using Eqs. (|2I3I14P in Eq.([n]) leads to: 

Fig. 5a shows the agreement of Eq. (p!5|) with the simula- 
tions. When Es < J, the critical size will be larger than 
4, leading to deviations from (p!5|) [ll]. 

In regime II, the thermodynamic theory of nucleation 
applies [12], and the formation rate of holes reads: 

J = Por+e(^j e-«=/^ (16) 
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FIG. 5: tsat (o), and ti/2 (□) -the time at which reaches 
1/2- as a function of Es/ J and J/T, with ICl. (a) Low tem- 
peratures, varying J/T with Es/J = 1, and varying Es/J^ 
with J/T = 5. The sohd hues correspond to Eq.([T5]) with a 
prefactor of 1.5 (b) High temperatures, varying Es/J, with 
J/T = 2.5. The sohd hne corresponds to Eq.([T8]) with a pref- 
actor of 25, and 7e//, while the dashed hne corresponds to 
Eg.dTSl) with 7. 



where Gc is the value of G at the critical island size Sc, 
and G'^ — 9ssG|s=s^, with s = 7ri?^, F+c is the rate 
at which atoms attach to the hole when R = Rc^ and 
po ^ 1/a^ is the monolayer density. The concentration 
in the vicinity of the holes at the critical size should be 
at equilibrium, i.e. c = ciq ~ e~^^/^, so that F+c = 
2i\RcCe}(D/a). Thus, Eq.([T6]) leads to: 

J = ^ £;l/2^-l/2^(-2J-7r7'aV£^5)/T^ ^yj^ 

Upon substitution in Eq (pT]) , and using (|4|6p we find: 



^-4/3^2/32^1/6 

^^/"^(e^^/^-l)^ 



,(8J+7r72aV£^s)/3T 



(18) 



We have performed simulations at J/T = 2.5. From 
the Ising model [10], we find 7/J ^ 0.27. The line ten- 
sion 7 at such high temperatures must be corrected in 
order to account for the non-ideal character of the 2D 
gas of adatoms around the cluster. Using the results of 
Ref.[i3, we find: 7e///J ~ 0.42 0. The high tem- 
perature limit is obtained when i?c ^ 1, i.e. when 
Es/J < n^^'^a^^eff / J ^ 0.74. The simulation results 
with 0.12 < Es/J < 1 reported on Fig |5] indeed agree 
with Eq.([l8l) for small Es. 

Orders of magnitude may be obtained for the case of 
Au/graphite. Using J ^ 0.5eV[U], Es - O.lel/^, 
and a diffusion barrier ~ 0.24eF[16], we find that we 
are in regime II, with tsat ^ 10 min, and A ~ 5 x lO^a 
at 800 JC. Moreover, the de wetting of a nanometer- thick 
Ag/Si(lll) [8], lead to patterns similar to those obtained 



here in regime II, suggesting that our finding may be rel- 
evant for films with a small number of atomic layers. 
Quantitative experimental analysis of the dewetting pro- 
cess in various systems is still lacking though. 

We have focused on the formation of the bilayer islands 
because this process determines the break-up time tsat, 
and the initial lengthscale A. The subsequent dynamics 
combines processes which have been separately studied in 
the literature: layer-by-layer thickening of the islands [l7|. 
coupled to their sintering [18] and Ostwald ripening. We 
hope to report along these lines in the future. 
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At low temperatures, the equilibrium island shape has 
a height heq a square base of lateral size Leg, with 

heq/ Leq — Es/ J ■ 

The equilibrium concentration with islands on top of the 
-2J/T ^ Thus, the attachment rate 



monolayer is Ceq 



of atoms to kinks is ^ Dc^eq <^ r2, the kink 

atom detachment rate. Therefore, backward kink motion 
is negligible. 

[21] The results of Ref.[13] are obtained in a strictly 2D 
model. The thermal roughness of the layer, which is rein- 
forced when Es / 0, is neglected here. Following Ref.[13j], 
one essentially needs to add a multiplicative factor to 7. 
Therefore, this should not strongly affect the tempera- 
ture dependence of Eq.(|6|). 



